Dissipationless behavior of asymptotic non-Markovian dynamics within structured 

environments 
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We investigate the long-lived quantum coherence in the non-Markovian process of spontaneous 
decay of a two-qubit system subject to a local environment. The asymptotic population of the 
qubit is linked to the spectral density of the reservoir through a general functional relation between 
them. The dissipationless behavior of asymptotic processes, including the steady coherent dynam- 
ics with self-sustained oscillations and ultimate persistence of quantum correlations, are elaborated 
with respect to spectral parameters of the specified reservoirs, e.g., the Ohmic class and the pho- 
tonic crystal. We expect these results to contribute towards reservoir engineering with the aim of 
enhancing stationary quantum coherence in noisy environments. 
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Understanding the non-Markovian dynamics of dissi- 
pative processes induced by memory effects of the en- 
vironment is a fundamental subject of open quantum 
systems [H, 0]- In general, the non-Markovian process 
differs from the Markovian one by exhibiting the non- 
exponential decay or even dissipationless behavior; it be- 
comes evident in cases of low-temperature environments 
and strong system-reservoir couplings and the resulted 
dynamics of the system depends dramatically on the 
spectrum structure of the environment. The relevant 
physics of the non-Markovian dissipative dynamics in the 
strong-coupling regime has been studied in the field for 
many years (3|-|5|- Due to the complexity of multiple en- 
ergy scales of the overall system, e.g., that of the system, 
the environment, and of mutual coupling between the 
system and the environment, it is generally a challenging 
task to characterize the relevance of the system dynamics 
with the reservoir spectra. 

Intensive interests upon the non-Markovian dynamics 
have been stimulated by the progress of quantum in- 
formation science [6] since the task of decoherence con- 
trol, aiming at performing quantum manipulation with 
high accuracy under noisy environments, requires that 
the memory effect of the environment needs to be taken 
into account. The back action of the non-Markovian en- 
vironment will prolong quantum coherence and even lead 
to distinct longtime behavior for the system from Marko- 
vian processes. For instance, partial inhibition of radia- 
tive decay of atomic excited states can occur in struc- 
tured environments of photonic crystals |7','8'| and high-Q 
cavities [S], and the associated potentiality of entangle- 
ment preservation under noisy environments has been un- 
veiled in recent literatures ^M,- Modifying the property 
of the reservoir to reach the non-Markovian regime was 
shown to be practicable in many physical systems, e.g., 
in photonic crystal materials [111 ] and in hybrid systems 
composed of quantum dots or trapped ions coupled with 



optically confined ultracold atomic gases [T^] • It is hence 
a fascinating task to describe the character of the long- 
time dynamics in connection with the spectrum struc- 
ture of the non-Markovian environment, for the purpose 
of reservoir engineering to enhance stationary quantum 
coherence in noisy environments. 

In this Letter we expose explicitly the link between 
the asymptotic dynamics of the system and the spec- 
trum structure of the non-Markovian environment for 
the spontaneous emission model. We establish a general 
functional relationship between the asymptotic excited- 
state population and the spectral density of the reservoir. 
This enables us to carry on detailed analyses upon how 
to engineer the spectral parameters optimally to achieve 
stationary coherent dynamics of asymptotic processes for 
various reservoirs, e.g., the Ohmic class and the photonic 
crystal. Furthermore, we probe the persistence of quan- 
tum correlations-by notions of entanglement and quan- 
tum discord-of a two-qubit system undergoing local non- 
Markovian dissipation. The trapping phenomena of the 
two quantities and their dependence on the spectral den- 
sities of the reservoir are displayed. 

The model we are considering is composed of three 
parts: qubits A, B, and a local environment E coupling 
with A. Initially the environment is in the vacuum state 
and the total system is described as ptot(O) = Pab{0) <8i 
|{0fc})({0fc}|. In the rotating- wave approximation, the 
interaction between A and E assumes the form 

Hae = iOQ(T+a^+'^[oJkblbk + {gkbkcr++glblcr^)], (1) 

k 

where ujq is the transition frequency between upper and 
lower levels |±) of the binary system A and (j± are the 
raising and lowering operators. The field modes of the 
reservoir are described by the creation (annihilation) op- 
erators bj, (bk) with frequencies uJk and gk account for 
the coupling constants of them with A. It is known that 
the time evolution of the reduced density operator of the 
qubit A is given by [l| 
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where the amphtude c{t) satisfies the integro-differential 
equation 



c{t) + iujoc{t) + / c(r)/(t - T)dT = 0, 



(3) 



and the kernel function of the integral is given by f{t~ 
r) = J(c^)e-^"(*-^)da; with J{u;) = Efc \gk\^Siu;-ujk) 
representing the spectral density of the reservoir. For 
particular forms of J(w), e.g., the spectra of photonic 
crystals, the above equation allows a solution with non- 
vanishing c{t) in the longtime limit. Notwithstanding, 
resolution to the system of previous methods, either via 
numerical calculations [l^ or analytical approaches (say, 
by the Laplace transformation 0, [iJl and by replac- 
ing the reservoir modes with pseudomodes [isl Il6j). was 
mainly focused on the time evolution of the population. 
Straightforward revelation upon the connection between 
the asymptotic qubit population and the reservoir spec- 
trum has only been obtained for very few cases with par- 
ticular forms of the spectral densities d, . 

To this goal, let us consider the eigensolution of the 
Hamiltonian of Eq. ([T]) in the single-excitation sector: 
HaeI^bs) = E\^Bs), where the energy E is determined 
by the secular equation 
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and the eigenstate takes the form — b\+,{Ok}) 

X]/c bk\ — , Ifc) with the coefficients obtained as 
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Note that the described issue, ac- 



counting for the coupling between the discrete levels and 
the continuous spectrum, identifies a standard Friedrichs 
model [17] and the specified eigenvalue problem could be 
retrospected to the work by von Neumann and Wigner 
[iSi]. With the presumption that J{uj) > holds for 
w > 0, Eq. Q allows at most a single real root that 
is non-positive |19|] . For the mentioned environment of 
photonic crystals, the existence of the solution of atom- 
photon bound states was known well 7, 8] and the secular 
equation Q might allow multiple real roots in view that 
the medium could possess several band gaps. Here we 
consider the case that Eq. ^ possesses only a unique 
real root. In this direct connection between the 

asymptotic population and the reservoir spectrum can 
be established in view of the facts that: i) Eqs. 
and ^ determine a conclusive relation of b with J(w): 
b — 6[wo, J{uj)]; ii) the asymptotic population on the ex- 
cited state is given by: |c(t)p |cooP = b^. 

To make clear about the result (ii), we express the 
excited state as |+,{Ofc}) = b\^Bs) + ^I^'d(O)), where 
|$d(0)) represents the projective state of |+,{Ofc}) over 
the complemental subspace orthogonal to \^bs) and 



b = Vl — 6^ is the corresponding probability. The dy- 
namics generated by the Hamiltonian ([T]) yields that 
\^D{t)} = e^*'^-*^*|$£)(0)) which will decay entirely as 
i — >■ (X). So the amplitude of the excited state in the 
longtime limit is contributed solely by the ingredient of 
it in |$_Bs), which leads promptly to that |coo| = b^. Note 
that this fact has ever been revealed in the literature Q 
to describe incomplete decay of an atom in photonic band 
gap mediums. 

The above disclosed facts (i) and (ii) suggest an im- 
plicit functional relationship between the asymptotic 
population and the spectral density function, |cooP = 
b'^[ujo, J{u!)], which renders a peculiar perspective to look 
into the character of long-lived quantum coherence of the 
non-Markovian dynamics. It enables one to characterize 
the dissipationless behavior of pA and pab hi^ connec- 
tion with the details of any specified spectrum structure 
of the environment. Later on we will elaborate the is- 
sue by focusing on two sorts of typical reservoir spectra, 
the Ohmic class J{uj) = Jo{oj) and the photonic crystal 
J{uj) = Jp{uj) [see Eqs. ([5]) and dS]), respectively]. Before 
presenting these details, we highlight that the nonvan- 
ishing Coo suggests an attractive limit cycle [2d|, i.e., a 
close dynamical trajectory of pA with self-sustained oscil- 
lations as time approaches infinity. Because a stable limit 
cycle is definitely a phenomenon of coherent dynamics-it 
could never happen in the statistic equilibrium of Marko- 
vian processes-its peculiar existence can be viewed as a 
representative of the non-Markovianity of dissipative dy- 
namics. A visual description of coherent evolution of the 
limit cycle is shown in Fig. 1, where the trajectory of pA 
in Bloch space is illustrated for the dissipative process 
with an Ohmic reservoir. The orbit is shown spiralling 
into a stable horizontal cycle around the z-axis. The ra- 
dius of the limit cycle is determined by = |coo|a(0) 
where a(0) = 2\p^ (0) | is the module of the horizon com- 
ponents of the initial state vector of pa(0). 




FIG. 1: Trajectory of time evolution of pAit) in the 
Bloch space starting from an pure initial state |^/)a(0)) = 
cos(7r/12)|-|-) -I- i sin(7r/12)|— ). The spectral density of the 
reservoir is specified by J{uj) = Jo{oj) [see Eq. ([6])] with pa- 
rameters s = 5.5, uJc = O.Sloo and rjo = 0.08. 
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From Eq. ([4]), it is evident that E will descend mono- 
tonically as the dissipation strength [e.g., rj^.p in Eqs. 
([6]) and ^ responsible for the reservoir of Ohmic class 
or photonic crystals, respectively] increases. The alter- 
ation of b in Eq. ([5]) , relying on its contained integrand, 
is determined by the competition between the numera- 
tor and the denominator, both terms ascending with the 
dissipation strength. As a consequence, the asymptotic 
amplitude |coo| is generally not a simple monotone of the 
dissipation strength, but their relevance conditioned to 
details of J{uj). 

Let us consider the reservoir described by the widely 
used Ohmic class spectra 



(6) 



where Wc is the cutoff frequency and s is a parameter 
whose scope, s < 1, s — 1, s > 1, corresponds to sub- 
Ohmic reservoirs, Ohmic and super-Ohmic reservoirs, re- 
spectively. The secular equation Q in this case allows 
an eigensolution of the bound state as long as the param- 
eters satisfy r]~^ < ^r(s) [2l|, where T{s) is the gamma 
function. The functional relation of Eqs. (|H) and ([5]) 
gives rise to 



b(r]o,s,u>Q/uJc) 



1 



r]oX°e~ 
{x — kY 



-1/2 



(7) 



where k = E/uc is determined by the transcendental 
equation uiq/lOc — f!^ n^^Ll — — Numerical calcu- 
lation to the latter equation is needed to acquire k for 
specified parameters [rjo, SjUJq/uJc), with which the exact 



population 



— b can be derived from Eq. 



We outline in Fig. 2 the characteristic of the popula- 
tion I Coop with respect to different zones of the spectral 
parameters. In the case of low Wc/wg, a high value of 
'r]oT{s) indicates that the scope of the parameters (ryo,s) 
is relatively narrow to achieve a nonvanishing c^c- Note 
that there is a physical constraint of the coupling strength 
in order to validate the rotating-wave approximation for 
the model Hamiltonian ([Ij (reasonably rjo < 0.1 owing 
to \gk\ ^ wq). The dependence of jcoop on rjo is shown 
in Fig. 2(a) with lOc/ujo = 0.3. Our calculation displays 
that high Ohmicity about s > 5.25 is required, which 
may challenge the technology of the reservoir engineer- 
ing. For the situation with high cutoff frequencies, the 
dependence of the population jcooP on the parameter s is 
shown in Fig. 2(b). Note that in the limit of ojo/uJc ~> 0, 
the transcendental equation defines an implicit function 
K — nijio, s), hence b = b{rio, s) according to Eq. ([7]). For 
•qo = 0.08, the maximum of the asymptotic population 
jcooP — 0.9 is achieved at s ~ 2.34. 

For the second example we consider the structured en- 
vironment of photonic band gap mediums with the spec- 
tral density given by [1, 



Jp(a;) 



3/2 



= e(w - We), 



(8) 



0.35 



Poo 0.31 



0.27 




1.0 



Po. 0.6 



(b) 




FIG. 2: Long-lived population Poo = |ccx)P in relation to the 
parameters of Ohmic class spectra, (a) Poo as a function of 
r]o with Uc = 0.3a;o, s = 5, 5.25, and 5.5, respectively. The 
green-dotted line figures out the boundary determined by the 
critical condition rjoT^s) — ujo/oJc- For s = 5.5, the maximal 
Poo = 0.33 is achieved at -qo = 0.08. (b) Poo as a function 
of s with rjo = 0.08, ujc/cjo — 15, 50 and oo. In the limit 
cjc/i^o — >■ oo, the maximal Poo = 0.90 is achieved at s = 2.34. 



where is the band edge frequency, O is the Heavi- 
side step function, and rip accounts for the dissipation 
strength We have neglected the influence of the 

lower band by considering an atom with the transition 
frequency cuq around U!e- The specified integrals in Eqs. 
dH) and ([5]) for this case can be worked out analyti- 
cally. In brief, Eq. (|4]) reduces to an algebraic equa- 
tion {ujQ — E)y^u}e — E = -qpLoV^ and the corresponding 
real root determines the bound-state energy E = E{r]p). 
Consequently, according to Eq. ([5]) one obtains 



3/2 



2{lo, - E f/^ 



E 



2{uoe - E) 



(9) 

where the identity of the i?(?7p)-related algebraic equation 
has been used to achieve the latter equality. For fixed 
and the peculiarity of the relevance |cooP = b'^{rjp) is 
then clearly elucidated: i) for the case that the atomic 
level is inside the reservoir continuum (wq > We), is an 
increasing monotone of rjp; ii) for the atomic level outside 
the continuum (wq < We), & is a decreasing monotone 
of rjp] III) for the exact resonance case with cjq = We, 
the population becomes a constant |cooP = 4/9, which 
recovers the result of the previous literature Q . 

Let us consider further the dynamics of the two-qubit 
density operator pab and the ultimate persistence of 
quantum correlations kept in the system. Straightfor- 
wardly, the time evolution of p^B under the local dissi- 
pation can be described by the Kraus representation 



PAB{t) = ^ri(t) ® IbPAB (O)rl(t)0/B, (10) 



where 



ri(t) 







c{t) 
1 



(11) 
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The validity of the above result is justified by the fact 
that the expression of PA(t) in Eq. 1^ can be recast as 
PAit) = T,'i=l'^^{t)pA{0)Tlit). One may notice that c{t) 
will suffer a phase uncertainty for the longtime dynamics 
of the asymptotic process. It happens that the phase 
factor of c{t) does not affect any measure of quantum 
correlations between A and B. This can be seen since if 
we replace c{t) by its module \c{t)\ in r2{t), the yielded 
two-qubit state will differ from the original pAsit) only 
by a local phase shift ua ® Ib with ua =(o'^ °) and 
if = aTg[c{t)]. 

We employ two different notions of quantum corre- 
lations, entanglement and quantum discord, to illus- 
trate the persistence of quantum coherence in the de- 
scribed dissipative process. Reported in Fig. 3 are the 
two quantities yielded from an initial Bell state sub- 
ject to a local reservoir of the photonic crystal. The 
amount of entanglement of the two-qubit state is fig- 
ured out via the measure of concurrence j23| : Cab = 
max{0, "\/A]"~X) j=2 v^l' where Ai are eigenvalues of the 
matrix pABO'y ® '^vPab'^v '^y in descending order. The 
discord of pab, according to definition of Ref. [13, is given 

by QaB = TailHAkyJ2kPkS{p%\{Ak}) + S{pA)-S{pAB), 

where S{p) = — tr(plog2 p) is the von Neumann entropy. 
The first term in the definition, so called as conditional 
entropy, involves an optimization procedure over the pro- 
jective measurement {Ak} and a tractable approach is 
yet absent to deal with it for general two-qubit states. 
Favorably, we show in the below that for any pure input 
states, the discord as well as the concurrence of yielded 
states can be achieved analytically. 



(a) (b) 
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FIG. 3: Asymptotic behavior of quantum correlations of an 

initial state |'i/'AB(0)) = -^(| H ) + | h)) subject to a 

local structured reservoir of photonic crystals, (a) Quantum 
discord Q'ab a-s a function of r]p with Ue/i^o = 0.8, 1.0, and 
1.2. A constant Qab ~ 0.56 is obtained as oje ~ ^j^o- (b) 
Time evolution of the discord as well as the concurrence for 
the case of uje ~ i^o, where rjo = 0.08 is assumed. 

Without loss of generality, we express the pure initial 

state as 

IV'ab(O)) =a|-l-,^+)+/3|-,V5-), (12) 



where a and (3 are positive real numbers satisfying + 
= 1; |</?±) are two normalized states of the qubit B 
but not necessarily to be orthogonal. According to Eq. 
PTI| . the state will evolve as 

PAB{t) ^ a^C^{t)\~,^+){-,^+\ + \(j>AB){<l>AB\, (13) 

where c{t) = \/l - |c(t)p and I^ab) = ctc{t)\+, ip+) + 
— , (/?_). The concurrence is then obtained as CAB{t) = 
\c{t)\CAB{Q) with Cab{Q) = 2aP\{^-\ay\^X)\ denoting 
the concurrence of |V'ab(0)). By noticing that the state 
PAB{t) is of rank two, we can work out its discord through 
a scheme proposed in Ref. l25|. The result reads 

QAB{t) - h{\) + h{\A) - HXab) = QAB{c{t)), (14) 

where h{x) = — xlog2 x — [1 — x) log2(l — x) is the bi- 
nary entropy function. All the contained parameters. A, 
Ayi, Aab, hence the discord itself, are functions of the 
amplitude c. Their detailed expressions are shown as: 
A(c) = i[l + ^l-<?C\sm, Xab{c) = Aa(c), and 



Aa(c) = 2 + V ^2 ~ + \al3c{^+\^^W. (15) 

In the limit i — > cx), a steady |coo| is reached, hence 
the discord QAB{t) Qab — Qab['^('^)] (^^ ^^^^^ 
concurrence C^). For the specified spectra of photonic 
crystals, the relevance of Qa'b with respect to rjo as well 
as the time evolution dynamics are manifested in Fig. 
3. We point out that both the discord and concurrence 
of the output state (|13|) are monotonic functions of the 
upper-level population |cp, thus their longtime behavior 
is similar to that of |cooP described previously in connec- 
tion with the parameters of various reservoir spectra. 

In summary, the asymptotic non-Markovian dynam- 
ics of the spontaneous emission model was investigated. 
The dissipationless behavior of the model within various 
structured reservoirs has been elaborated via the func- 
tional relationship between the asymptotic excited-state 
population and the spectral function of the reservoir. 
Moreover, we have figured out the long-lived quantum 
correlations of a two-qubit system subject locally to the 
specified dissipative process. We expect these studies to 
contribute towards engineering of the reservoir spectra 
that aims at attaining steady quantum correlations un- 
der noisy environments. 
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